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Abstract 

We consider the motion of tachyons (faster-than-light particles) 
in the framework of General Relativity. An important feature is the 
large contribution of low energy tachyons to the energy-momentum 
tensor. We also calculate the gravitational field produced by tachyons 
in particular geometric arrangements; and it appears that there could 
be self-cohering bundles of such matter. This leads us to suggest 
that such theoretical ideas might be relevant to major problems (dark 
matter and dark energy) in current cosmological models. 
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1 Introduction 



While there is no reliable evidence for the existence of faster-than-light parti- 
cles (tachyons) some physicists have considered them a theoretical possibility. 
[1, 2] Yet it appears that many carry a prejudice against the legitimacy of 
tachyonic theories because of concerns about causality. This comes from a 
familiar "paradox" that posits communication via tachyons that appear to 
travel backwards in time. In Appendix A I have tried to dampen that preju- 
dice by looking carefully at the assumptions behind the paradox. If tachyons 
interact only weakly with ordinary matter, yet are profusely present through- 
out the universe, then it may be that those objections become nugatory. 

The present paper undertakes an investigation of how faster-than-light 
particles might be considered within the General Theory of Relativity. There 
has been some earlier study of this [3] but the present work covers consid- 
erable new ground and perhaps may make some interesting contributions to 
current models of cosmology. 

We start, in Section 2, with how they would be described in Special 
Relativity, noting the particular role of the energy- momentum tensor T^ u . 
Then we turn to General Relativity. In Sections 3 and 4 we study the motion 
of a tachyon in two familiar gravitational models. Then in Sections 5 and 
6 we derive the gravitational field produced by a uniform flow of tachyons; 
and find that there is an attractive force on other nearby tachyons which 
increases in strength as their energy decreases. In Sections 7, 8, 9 and 10 
we explore various models for gravitational fields produced by static flows of 
tachyons, using the linear approximation to Einstein's equation; and we find 
that self-cohering bundles of tachyonic matter are feasible. Section 11 has 
a note about how tachyons might contribute to gravitational lensing. These 
results, summarized in Section 12, suggest that tachyons should be given 
closer attention in the ongoing study of cosmology. 

2 Special Relativity 

In Special Relativity, we have the motion of a particle described by the four- 
vector coordinate £ m (t), fi = 0, 1, 2, 3 and r is a parameter called the proper 
time. For a free particle moving in a straight line with velocity v we have, 

e(r) = ( 7 r, 7 vr). (2.1) 
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Then, writing £ M = d^/dr, we calculate 

V^r^"^7 2 (l-A (2.2) 

where we use the metric of Special Relativity, rj^ (with i] q = +1 and the 
other diagonal components equal to -1) to raise and lower indices. We define 
7 = 1/ \J\1 — v 2 \ so that this is always a real number. Then we see that for 
normal (slower-than-light) particles the invariant in Eq. (2.2) has the value 
+ 1 and for tachyons it is -1. (I am using units in which the velocity of light 
is c = 1.) 

Trying to follow conventional habits we can define the "energy" and "mo- 
mentum" of a free tachyon as follows: 

E = m 7, p = m 7 v, E 2 — p 2 — m 2 (2.3) 

where m is a (real, positive) mass parameter. This allows us to characterize 
the full set of free-particle tachyon states as, 

E > 0, p 2 > m 2 ; (2.4) 

and the direction of the vector p covers the entire sphere. 

Troubles occur when one regards this E and p as a 4- vector. For tachyons 
this is a space-like vector and so the sign of the first component can change 
under a proper Lorentz transformation. This can lead to unnecessary confu- 
sion about negative energy states. This confusion was straightened out in an 
earlier study [4], where it was shown how to properly use the labels for "in" 
and "out" states in any interaction. It was also noted there that the best way 
to talk about energy and momentum is by starting with the familiar second 
rank tensor. 

The energy-momentum tensor for a classical particle can be written as 

T^(x) =mj dr 5\x - f (r)) (2.5) 

and we note that the energy density, T 00 , is always positive. For the free- 
particle motion noted above, we calculate that this tensor density has the 
diagonal elements p, Pu, P22, -P33, where 

p = m7; p.. = rwyv 2 . (2.6) 
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Note that for tachyons the pressure terms (X), Pa) come out larger than 
the energy term (p); this is exactly the opposite of the situation for normal 
particles. 1 

To illuminate this unfamiliar situation, imagine the neutrino as being 
a tachyon. Then its contribution to the energy-momentum tensor, in the 
pressure term, would be larger than what one would ascribe to a massless 
particle by a factor of v 2 = (p/E) 2 = (E 2 + m 2 )/E 2 . If one were to assume 
a mass of around 0.1 electron Volts and an average energy corresponding to 
the cosmic background temperature of around 3 °K, then this would give 
a numerical factor of v 2 ~ 10 5 . That would, presumably, have a significant 
impact on current calculations of cosmological models. [5] 

For some further musings on kinematical aspects of tachyons, see Ap- 
pendix B. 

3 General Relativity metric for uniform cos- 
mology 

The Robertson- Walker metric, describing a homogeneous isotropic universe, 
is 

±ds 2 = dt 2 - A(t)[B(r)dr 2 + r 2 d6 2 + r 2 sin 2 6 dip 2 ]. (3.1) 

where we have identified the four (spherical polar) coordinates t, r, 6, ip. The 
functions A(t),B(r) are commonly written as: 

A{t) = a(t) 2 /c 2 ; B(r) = 1 _^ 2/jR2 - (3.2) 

With this metric we construct the following equations for the geodesic. 

i+ -A'[f 2 B + 9 2 r 2 + <p 2 r 2 sin 2 d] = 0, (3.3) 

Z 

-ABr + r[-iA'B - rAB'] + \A\r 2 B' + 9 2 2r + ip 2 2r sin 2 9] = 0, (3.4) 

-Ar 2 6 + 6[-iA'r 2 - r2Ar] + ]-A[p 2 r 2 2 sinO cosO] = 0, (3.5) 

z 

-Ar 2 sin 2 e (p + ip[-iA'r 2 sin 2 6 - rA2r sin 2 9 - 9Ar 2 2 sinQ cosO] = 0, (3.6) 



1 It is usually so that a repeated index implies summation; but there are exceptions, as 
seen here. 
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where the "prime" on A or B means derivative with respect to its argument. 
We can integrate the last two equations, yielding 

Ar 2 sin 2 6 tp = M, (3.7) 
A 2 r A 6 2 + M 2 /sin 2 6 = L 2 , (3.8) 

where M and L are some constants. Then, substituting these results back 
into the first two equations we also get the integrals: 

A 2 Br 2 + L 2 /r 2 = 2£, (3.9) 
i 2 -2£/A = Q, (3.10) 

with two more constants, S and Q. 

Using these results we now calculate the invariant 

g^x»x» = Q, (3.11) 

So, for normal (slower-than-light) particles we would choose Q = +1; for 
light we would choose Q = 0; and for tachyons, we choose Q — — 1. We 
shall see that in some instances the constant 28 may be taken as v 2 for some 
initial velocity v, while L and M relate to the angular momentum for unit 
mass. 

From Eq. (3.9) we see that motion of a tachyon will be contained in 
the case of a closed universe, (e = +1), even though its velocity may be 
arbitrarily large. This follows because the factor B(r) becomes infinite at 
r = R, and thus the radial coordinate r must stop its increase at that point. 

Another simple solution of Eq. (3.9) is a circular orbit (r = 0) at a radius 
r = L/v. But such orbits are not stable: for B > they will spiral outward; 
and they will spiral inward for B < 0, which may occur for r > R in a closed 
universe. 

4 Tachyon in Schwarzschild metric 

Outside of a central source of gravitation, we have the Schwarzschild metric, 
±ds 2 = A(r)dt 2 - [A(r)-^r 2 + r 2 d6 2 + r 2 sm 2 9 d V 2 }. (4.1) 
where A(r) = (1 — r s /r) and r s = 2GM. 
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The solution of the geodesic equations in this metric yields the following. 

rV = £, At = -y, A-\r 2 - 1 2 ) + L 2 /r 2 = -Q (4.2) 



where 7 and L are constants of integration and Q is the same constant (with 
values +1, 0, -1) defined in Eq. (3.11). 

For slower-than-light particles (Q=+l) this gives the Kepler orbits for 
r > r s , along with small relativistic corrections. For tachyons (Q=-l) we see 
that this does not give localized orbits but only scattering states. 

Introducing the variable, u = 1/r (</?), and introducing the asymptotic 
velocity v, with 7 2 = Q + v 2r y 2 , L = bwy, we get the equation for the orbital, 



From this equation, we can calculate the scattering angle to first order in 
the gravitational field strength and find the result 



For light (Q=0), we have x — 2, which is the well known result. For tachyons 
(Q=-l), x nas the value 2 at high energies and then decreases to 1 as the 
energy drops. For ordinary particles (Q=+l) at very high energies (v close 
to 1) x is again 2; and at much lower velocities we should not accept the first 
order approximation but solve the exact Kepler problem. 

We can also easily look at a "head-on" collision by setting L = in Eq. 
(4.2). Then we see how a tachyon will behave approaching the edge of a 
black hole, r — >■ r s . 

5 Metric in cylindrical symmetry 

For tachyons as a gravitational source, the simplest model is a straight 
line of particles, flowing constantly along the z-axis, with equal velocities 
in both directions. With cylindrical symmetry we use the four coordinates 
£ M = (t,r,0,z), where r = yjx 1 + y' 2 . The source is then a diagonal energy- 
momentum tensor density, following Eq. (2.6), 



T 00 = P1 5(r)/(2nr), T n = T 22 = 0, T 33 = pjv 2 5(r)/ (2nr), (5.1) 



u' 2 + u 2 - 



Qr s u + w 2 7 2 




(4.3) 



v Ar i z 



(4.4) 
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where p is the mass per unit length. We now make this ansatz for the metric: 

±ds 2 = A(r)dt 2 - B{r)dr 2 - r 2 d6 2 - C(r)dz 2 . (5.2) 
Next, we calculate the elements of the Christoffel symbol. 

(5.3) 
(5.4) 
(5.5) 
(5.6) 

where the prime means derivative with respect to r. 

From this, we calculate the components of the Einstein tensor G^ u = 

(5.8) 
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r 2 A" A' 2 

^22 — 2g l ~^ o A2 o a d ~ ~^ nrvz ~ o Art o dT^ J 1 ' ' " 

G »=^ { -T + ^ + £^-St + ^ } - (510) 

If we set all these components equal to zero, we find the solutions 

A = ar a , B = hr fi C = c r 7 (5.11) 

0=-^> (5.12) 
a + 2' a + 2 v 1 

where a is some constant (7^ —2). 

For small a we have the ("weak field") solutions: A = 1 + a ln(r), B = 1, 
C = 1 — a Zn(r). 

Now, using the Einstein field equation = —8nGT^ u , with Eq. (5.1) 
we find a 12 = SnGp^ / (2n). This leads to the result 

Force = -T^ = -AGpj/r (5.13) 
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for the case of a slow particle being attracted by a line of mass-density p. 
This is to be compared with the result of Newtonian gravity; and there is a 
difference of a factor of 2. (This problem in General Relativity was studied 
long ago by Levi-Civita [6] and he obtained the correct Newtonian force.) I 
have no explanation of this discrepancy. Maybe it has to do with the line of 
matter extending to infinity. 

Similarly, for a fast particle near a line of streaming tachyons we find the 
attractive force 

Force = = -AGfrfv 2 (z 2 )/r. (5.14) 

What is novel here is that the strength of the gravitational attraction between 
tachyons increases as the velocity increases (without limit). 

There is some question about the legitimacy of this result due to the non- 
linearity of the equations and the singularity at r = 0. Please see Appendix 
C for resolution of this question. 



6 Circular model 

To make a more realistic model of tachyons as a gravitational source, we can 
imagine something like a circular flow, which is limited in its spatial extent. 
The full solution of Einstein's field equations then becomes more difficult; 
but we can start out by making the weak field approximation. The general 
approach (see, for example, [7]) goes like this: 

d 2 

(g- t - V 2 ) V = -IGxGT^, (6.1) 

9tiu = V^u + V _ h = rf v h^, (6.2) 

with the side condition (gauge choice) 

= °- ( 6 - 3 ) 

This gives us Newtonian gravity in the familiar case of mass at rest (Too = 
p(x)); and we use this now for a tachyon source, following (2.5), (2.6). 

If this is a large uniform and static circulating flow of tachyons and we 
look close to the edge of it, then we can approximate it as an infinite cylinder, 
along the local z-axis, and immediately see the result. From h 33 we get an 
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attractive force on a unit mass at rest; and it is without that previous extra 
factor of 2. 

F = -2GP/r (6.4) 

where r is the distance from the center line of this cylinder and P is the 
"pressure" (m'-fv 2 ) per unit length of the source. For large values of v, which 
means low energy for the tachyons in the source, this force becomes large, 
proportional to \v | (or l/E). This suggests remarkable new physics, if such 
things actually exist. 

Let us see how this force works upon another tachyon travelling nearby. 
The geodesic equation gives us 

where this 7 describes the free-moving tachyon. So, a low energy tachyon is 
accelerated toward the flow of other low energy tachyons. This effect increases 
with the velocity of each component and suggests how such a structure - a 
self sustaining and accreting mass flow of tachyons held together by gravity 
- might evolve. 

For example, if we have a total mass M of tachyons moving within a tube 
of radius r at large velocity V in a circle of radius R, then an additional 
tachyon moving with large velocity along the edge of this tube will have its 
trajectory bent to keep up with that circle provided that GMV/r > ir. 

In Section 10 we shall carry this study further to include tachyon orbits 
spiraling around the center line of this source. 

Any such structure would also have gravitational effects upon ordinary 
matter and upon light beams. This is what the geometric theory of gravity 
requires and we now try to consider this in some generality. 

7 Static gravitational fields 

Now we want to follow the first-order approach of Eqs. (6.1 - 6.3) for a general 
distribution of ordinary matter and tachyons that is static, not changing in 
time. We would like to be able to specify a source T MJ ,(x), for example, one 
that has a simple structure like that of Eq. (2.6), where the energy density 
p and the pressure terms Pa are now limited by some spatial envelope we 
might choose. 
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This model has an obvious shortcoming in that d^T^ is not zero. (It 
would be zero if we only had the energy density and ignored the pressure 
terms; but for tachyons the motion is the main thing.) At first, I thought 
to put that criticism aside with the following argument. The particles that 
comprise this source experience some acceleration in order to remain in a 
confined region and we assume that the requisite forces are themselves grav- 
itational. Thus the error introduced into the Einstein equation will be of 
second order in the constant G and we shall ignore it. (This seems analogous 
to the distinction between the ordinary derivative of and its covariant 
derivative.) 

This argument, however, is false on two fronts. In terms of the mathe- 
matics: The Einstein equation = SnGT^ has, by construction, a left 
hand side that automatically satisfies the conservation law. Thus, using a 
right hand side that fails this test may (will) lead to incorrect solutions. In 
terms of the physics: For non-relativistic particles moving in gravitational 
orbits we know that the average potential energy is just as big as the kinetic 
energy. Therefore, it would be wrong to ignore the former and only keep 
the latter in constructing the energy-momentum tensor. (Again, for non- 
relativistic particles these pressure terms are generally negligible, of order 
v 2 /c 2 compared to their rest-mass; but for tachyons they dominate.) 

First, we shall find a general statement about the gravitational fields 
produced (in first order) from a static source that is localized in space. Then 
we shall proceed to some particular models. 

We have both T^ u and independent of the time variable, with the 
former assumed confined in some region of space. We also take the time- 
space components of each of these tensors to be zero. The basic equations 
are, 




(7.1) 
(7.2) 

(7.3) 
(7.4) 



dihijtx.) = d i T ij (x) = 0, 



with i, j = 1, 2, 3. The asymptotic form of Eq. (7.3) is 




(7.5) 
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but from Eq. (7.4) we calculate, 








(7.6) 



This means that we get no long range force from the spatial parts of the 
source, as suggested by Eq. (7.5); that comes only from the T 00 component, 
Eq. (7.2). Is this well-known? For ordinary (slow) matter this is no big deal 
(maybe it gets involved when one asks about v 2 /c 2 corrections to Newtonian 
gravity). But for tachyons it would be the major feature of their gravitational 
role. If there are such static, localized configurations of tachyons as we have 
been speculating here, then they may make major contributions to the local 
gravitational field; but they produce little or no long range gravitational 
attraction. What a lovely way to explain both dark matter and dark energy! 

Incidentally, Eq. (7.6) may be recognized as a generalization of the Virial 
Theorem, which is well known in non-relativistic mechanics, either the clas- 
sical or the quantum variety. 

What can one say about the next term, after ~ 1/r in the asymptotic 
expansion of Eq.(7.3)? It would appear to be ~ 1/r 2 with coefficients 
/ d 3 x Xk Tijix). These can be shown to vanish, however, by considering 
the identities / d 3 x x k xi diT^ = 0. The leading term will then be of order 
1/r 3 , which is what one expects of a quadrupole source. 

In what follows we shall try to construct models of the energy-momentum 
tensor that obey the conservation law, are static in time, and localized in 
space. That means they need to involve two parts which we may call the 
"kinetic energy" part and the "potential energy" part. 

8 Circular flows 

We shall be interested in problems with 3-dimensional and also 2-dimensional 
symmetry. So let us do this analysis first for n-dimensional Euclidean tensors 
Tjj(x), with i,j = 1, . . .n. We have the radial coordinate r = y/J2i 

For circular motion we can see that the kinetic part of the nxn matrix Ty 
should look like this 



1 (kin) 



(5ij - XiXj/r 2 )b(r) ~ (vi Vj) 



(8.1) 



where any multiplying envelope b(r) should be non-negative. 
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The easiest way to verify that this does represent circular flow is to cal- 
culate XiT^ m ^ and see that it is zero, since there should be no velocity in the 
radial direction. 

It is then convenient to introduce the nxn coordinate matrices / and D 
as follows . 

lij &ij 5 ^ij ^ij TlXiXjjv , (8.2) 

Trace(I) = n, Trace(D) = 0, D 2 = (n - 1)1 - (n - 2)D. (8.3) 

Then, we construct the full tensor as 

T tj = 5 i3 a{r) + D^r) = (8.4) 
<ya(r) - (n - l)6(r)] + nT\f n) . (8.5) 

and thus we identify [a — (n — 1)6] as the "potential energy" portion of T. 
From the conservation requirement 8(1^ = 0, we get the constraint 

r n a'(r) = (n-l)(r n b(r))', ia - (n - 1)6] = n(n - l)b(r)/r > 0. (8.6) 

Hi I 

Thus the part we identify as potential energy in the tensor is diagonal and 
increases with r. Assuming that both the functions a(r) and b(r) go to zero 
at large r, this means that the potential energy is negative - an "attractive 
well." From this we also learn that r n ~ l dr a(r) = . 

One more general result is the following, regarding the n-dimensional 
Lapacian operator. 

V 2 D g(r) = D[g" + (n - l)g'/r - 2ng/r 2 }. (8.7) 

Then, if we write 

h ij = S ij f(r) + D ij g(r) (8.8) 

r n f (r) = (n- \){r n g{r))\ (8.9) 

where we have required dihij = 0, then we get the solution 

g ( r ) = -^f- [r~ n / s n+1 ds bis) + r 2 / s- l ds bis)] (8.10) 
n + 2 jo Jr 

f(r) = -16nG(n - 1) r J™ s^ds b(s). (8.11) 
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There are alternative solutions in terms of the source function a(r); but they 
are determined by the relations stated above between / and g, on the one 
hand, and a and b on the other. 

From the fact that we have b(r) > 0, we can describe the two functions 
qualitatively as follows: g(r) is everywhere negative and goes to zero at both 
r — > and r — > oo; f(r) has everywhere a negative derivative and falls off to 
zero as r increases, thus it must be positive at all finite r. 

The asymptotic solution at large r is, as expected, 

1fi7r(7 r°° 
g{r)~--—-r- n s n+1 dsb(s), (8.12) 
[n + 2) Jo 

and the function /(r) falls off even faster. 

From the equations above one can derive a number of identities, for ex- 
ample, 

poo <n 1 poo 

f(0) = -n(n-l) dr r- 1 g{r) = 16nG / drrb{r). (8.13) 

jo 2 Jo 

Finally, that portion of the metric which depends on the function g(r) 
can be simplified, as follows. 

n 

dxi dxjlSij - nXiXj/r 2 ] g(r) = (8.14) 

Y^(dxi) 2 g(r) - n (^Xi dx}j g(r)/r 2 = (8.15) 
(g(dxi) 2 - n(dr) 2> J g(r). (8.16) 

9 Spherical source 

We start by applying the machinery of the previous Section to the case of a 
spherically symmetric source, n — 3. We have the following solution of the 
field equations, with r = \/x' 2 + y' 2 + z' 2 , 

g (r) = -™KG[r- z / s 4 ds b(s) + r 2 / s^ds b(s)\ ~ / d 3 x r 2 b(r), 

5 Jo Jr 5r J 

(9.1) 

and the function /(r) falls off even faster at large r. 



12 



We still have 

h 00 = -AG J d 3 x'^^ = AV (r) ~ -r~\ (9.2) 

Next, we re-write the metric, changing from Cartesian to spherical polar 
coordinates. 

ds 2 = [1 + A]dt 2 + [-1 + B]dr 2 + [-1 + C] r 2 [d6 2 + sin 2 6 dip 2 ] (9.3) 
A(r) = 2V + ^f, B(r) = 2V -±f-2g, C(r) = 2V - \j + g. (9.4) 

Then, we calculate the Christoffel symbols and put them into the geodesic 
equations. We immediately integrate two of them, 

i = 7e" A , r 2 (p = Le c ', (9.5) 

where 7 and L are constants and, as usual, we choose 6 = tt/2. Now the 
radial equation, in strictly first order form, leads to 

iff + 7 2 A(r) + L 2 /r 2 = 28 = 7 2 - Q. (9.6) 

The critical question we want to answer is whether tachyons may find 
bound orbits in this configuration of circulating tachyons (plus stationary 
matter). The potential Vq is negative, falling off to zero at large r; and the 
additional component of A(r), which is f(r), merely adds to the familiar 
L 2 /r 2 term. Thus, with Q = — 1 for tachyons, we see that £ > and so 
conclude that this equation does not allow bound states. 

Since we are seeking a self-supporting arrangement of tachyons, held to- 
gether by their own gravitational forces, we must look for some other geom- 
etry. That is where we turn now. 

10 Cylindrical source 

We now apply the above machinery to the case of a cylindrically symmetric 
source, n = 2. Here we have T 3i = h 3i = 0; and we have the following solution 
of the field equations, with r = y/x 2 + y 2 : 

16 r r r°° 2G r 

g (r) = TxG\r~ 2 / s 3 ds b(s) + r 2 / s~ x ds b(s)] ~ / d 2 x r 2 b(r), 

4 Jo Jr r 2 J 

(10.1) 
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and the function /(r) falls off even faster at large r. 
We still have 

h 00 = -AG J rfV ^°°_^ = 4Vb(r) ~ ln(r/R), (10.2) 

/,33 = -AG J d 3 x'^^ = W 3 (r) ~ Zn(r/it>), (10.3) 

where i? is some measure of the linear extent of this source; in other words, 
these potentials are negative with a positive slope. 

Next, we want to calculate the Christoffel symbols and put them into 
the geodesic equation. Changing from Cartesian coordinates to cylindrical 
coordinates (t,r,6,z), we find, 

ds 2 = (1 + A)dt 2 + (-1 + B)dr 2 + (-1 + C)r 2 d6 2 + (-1 + D)dz 2 (10.4) 
A(r) = [2V + 2V 3 + f\, B(r) = [2V -2V 3 -g], (10.5) 
C(r) = [2V - 2V 3 + g], D(r) = [2V + 2V 3 - /]. (10.6) 

The non-zero Christoffel symbols, to first order, are then 



1 ir — 1 tt — 2 ' 


(10.7) 




(10.8) 




(10.9) 


r(l + S-C) + — , 


(10.10) 


r^ r = -r; 2 = - V 


(10.11) 



Three of the four geodesic equations are immediately integrated, as fol- 
lows, 

i = 7e" A , z = ne D , r 2 9 = Le c ', (10.12) 

where 7, k, L are constants, which might be named energy, linear momentum 
and angular momentum, respectively. 

The fourth equation, for the radial coordinate, is more complicated. In a 
strictly first-order approximation it can be integrated to 

(f) 2 + L 2 /r 2 + 7 2 A + k 2 D = 28 = 7 2 - k 2 - Q, (10.13) 
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where, as always, Q = — 1 for tachyons. 

The question we address now is whether this equation allows bound orbits, 
which would imply particles moving in spirals around the z-axis. The effective 
potential in this radial equation is, 

^ + U(r), U{r) = ( 7 2 + ^)(V (r) + V 3 (r)) + ( 7 2 - n 2 )f{r)/2. (10.14) 

For bound states we would require, at least, that U(r) be a function that 
increases with r at large r. Both potentials Vo and V3 are of this character; 
and the function /(r) falls off faster at large r. So this requirement is satisfied. 
(This is no trivial result. It depends on how various plus- and- minus signs 
have worked out in this calculation.) The next requirement for bound states 
is that the integration constant S have a value equal to or greater than the 
minimum point of the total effective potential. 

To get an idea about this, we can say that for a free tachyon, or one in a 
scattering state, we would have, 

7 2 = l/(t; 2 - 1), K 2 = i 2 v 2 z) (10.15) 

which leads to 28 = rf{v 2 — v 2 z ) > 0. 

However, the integration parameters 7 and k are free from such con- 
straints, and we see that S certainly can be negative, thus allowing for bound 
state solutions. Note the contrast between this result and that of the previous 
Section, which looked at spherical symmetry. 

This result expands on what we had seen in the earlier work of Sections 
5 and 6, allowing for the circulatory motion in the source, along with the 
translational motion, relative to the z-axis. We conclude that streams of 
tachyons could form into coherent localized bundles held together by their 
own gravitational forces. 

We are now in the position to invite speculation about detailed clustering 
of large numbers of tachyons (assuming that such things exist) in something 
akin to rope-like structures throughout the universe. How big or small might 
such ropes be - as tight bundles or wide tubes; how might they connect to 
themselves - in simple circles or in knots, in doughnuts or in halos; and how 
might they attach themselves, gravitationally, to large masses of ordinary 
matter? Model-making and detailed calculations are needed to work out the 
answers to these questions; and further expert analysis is needed to see how 
such models may comport with known observational data. 
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11 Gravitational lensing 



The deflection of light beams by the gravitational field of dense matter is an 
important tool in the observation of galactic structures; and it has been used 
to impute the distribution of dark matter. If tachyonic matter actually exists, 
and is collected in static bundles, as has been suggested in the analysis of this 
paper, then one would want to ask how that might contribute to gravitational 
lensing. 

In the standard analysis, one works from the assumption that the source 
of gravitational fields is slow moving massive matter; and so one starts with 
the Newtonian potential, which we have called Vo(x) and carries out a sort 
of ray-tracing calculation to get the pattern of light beam deflection. 

If tachyonic matter is a significant presence, then this calculation needs 
to be augmented. It turns out that one should use, 



where the 3-axis denotes the original direction of the light beam. For low 
energy tachyons, the T 33 could be quite large. 

12 Summary 

There are several new results noted in this paper about how faster-than-light 
particles (tachyons) would behave under the General Theory of Relativity. 

First is the form of the energy-momentum tensor T 7 ^, where the spatial 
components would be larger than the time component. This is in contrast 
with the situation for ordinary matter; and it raises the question of how this 
might effect theoretical models of the overall structure of the universe, as in 
the Robertson- Walker metric. 

Second is the attractive gravitational force between tachyons and the sug- 
gestion of a stable model for a closed stream of such particles. This leads 
to questions about how such structures, if they exist, might relate to astro- 
physical observations of significance in the context of current cosmological 
models. The possibility of tachyons explaining some or much of what is cur- 
rently called "dark matter" is an obvious thought, which I leave to better 
experts for consideration. 

Another significant result, noted in Section 7 , is the absence of a long 
range gravitational force (~ 1/r 2 ) produced by localized static bundles of low 




(11.1) 



16 



energy tachyons. This prediction is something that may be readily subjected 
to observational verifiction, to prove or disprove the substantial existence of 
such tachyons. A recent study [8] of tidal properties of satellite galaxies may 
be relevant. 

In Appendix B we have noted a few additional aspects of tachyon kine- 
matics that might lead to observational tests of whether or not these things 
actually exist. 

This paper has been about classical particles that move faster than light, 
in the context of General Relativity. Further theoretical study is needed 
to look into tachyonic fields, how they propagate and how they may be 
quantized. 
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Appendix A: Fear of tachyons 

Lets start by reciting a version of the most common paradox in consider- 
ing the possibility of faster-than-light particles (tachyons) within the Special 
Theory of Relativity. For further reading, see Chapter 9 of the review paper 
[1] and other references therein. 

A space Ship moves away from Earth (along the positive x-axis) with 
velocity v < 1. At time t A , when the Ship is at coordinate x , a tachyon is 
sent out from Earth with velocity V > 1 in the direction of the Ship. This is 
event A. That tachyon reaches the Ship (this is Event B) at a later time ts- 
This is the story in the Earth frame of reference; and it may be summarized 
with the following formula. 

t B -t A = Tp-. (A.l) 
V — v 

What is the story as seen in the space Ship frame of reference? Here 
we use the standard formulas for a Lorentz transformation and come to the 
formula 

t' B -t' A = 1 {l-vV){t B -t A ). (A.2) 

This says that, if vV > 1, then the tachyon arrives at the Ship before it was 
sent from Earth according to observers in the Ship frame of reference. 
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This is the first chapter of the paradox. It is no surprise to those of us 
with some sophistication with the Special Theory of Relativity. We know 
that if two events are separated by a space-like interval, then the actual time 
sequence of those two events may be different as seen in different Lorentz 
frames. 

The full horror of anti-causal happenings appears when one adds a second 
chapter to this story; observers in the Ship respond to receiving the tachyon 
by sending another tachyon back to Earth, even faster, and this is seen to 
arrive back on Earth before the first signal was sent out. 

In order to avoid jumping to the conclusion that tachyons cannot exist, 
we should ask what lesser restrictions upon them might be considered. The 
following is a list of assumptions used in the story recited above. 

a) Special Relativity is correct; 

b) Tachyons exist; 

c) Tachyons can be created and sent out at will; 

d) Tachyons can be detected reliably. 

It is assumptions (c) and (d) that we should address first, before we reject 
assumption (b). 

We can suggest that if tachyons exist, they should occur as particle- waves 
of the sort we describe in the quantum theory; they may be emitted and 
absorbed by interacting with ordinary matter but these interactions must 
be very weak. This could allow for a very uncertain experimental set-up 
regarding controlled emission and detection (absorption) of tachyons. One 
would want to make a quantitative theory and show that the probability 
of observing the sort of nasty "paradox" described above is swamped by 
unavoidable uncertainty or background noise. 

One simple example is the following. Suppose the emission of the tachyon, 
as described above, was governed by a wave equation and appeared to ob- 
server A as a wave packet such as this. 

(p(x, t) ~ e l{kx - ut) e - r(t - x/v) 6(t - x/V), (A.3) 

where V = k/u is the group velocity. This packet has a size, in the original 
reference frame, of Ax = V/T. If we perform the Lorentz transformation 
on these coordinates to see what this wave packet looks like in the frame of 
observer B, then we find that the size of the wave packet is Ax' = V/[jT(l — 
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vV)]. Thus, when we approach that critical condition when vV ~ 1, the wave 
packet is so large that it encompasses both the "emitter" and the "receiver" 
at the same time. This gives us a view of how wave properties for tachyons 
can indeed erase the distinction between "emission" and "absorption" of a 
particle that travels faster than light. There is no paradox here. 

One common idea is that neutrinos might be a candidate for tachyonity: 
they have weak interactions and are nevertheless plentiful throughout the uni- 
verse. Against this idea, one might argue that neutrinos cannot be restricted 
by disallowing assumption (c) and/or (d), above, because of the possibility 
(in fact an established experiment) of using an accelerator to create a beam 
of neutrinos, which can be pulsed at the emitter and then detected reliably 
at a distant absorber. To refute this objection, one should notice a particular 
condition for the "paradoxical" situations discussed earlier. In order to get 
the effect of time-sequence reversal, between the emission and the absorption 
of a tachyon, one must have the source and the receiver moving away from 
each other. If they are moving toward each other, then there is no such ef- 
fect predicted by the Lorentz transformation. (Notice what happens if you 
change the sign of v in Eq. (A. 2).) The arrangement of an accelerator send- 
ing out a large pulse of particles in the forward direction of the high energy 
beam involves just this second situation: source and receiver approaching 
one another; so there is no paradox to be expected here. 

Of course, other arguments may be advanced and need to be studied. 
But, at least at present, I believe this is a constructive step in clarifying this 
controversial topic. Tachyons may exist and should not be dismissed from 
consideration in either theoretical or experimental studies. 

Appendix B: Kinematical musings 

If we calculate the number of free tachyon states that are "on the energy 
shell" inside a spatial volume V, we get 

]T 5(E-E ) = ^E oPo (B.l) 

states 

which is exactly the same formula one has for normal (slower-than-light) 
particles. 

So, if one wanted to ask whether neutrinos might actually be tachyons, 
then the search for neutrino mass, as carried out by looking at the end point 
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of electron energy spectrum in beta-decay, takes on new possibilities. The 
phase space factor (B.l) then gives us the spectrum for electron energies E e : 



(E — E e ) for massless neutrinos (B.2) 



(E — E e )J (E — E e ) 2 — m 2 for ordinary massive neutrinos (B.3) 



(E — E e )y (E — E e ) 2 + m 2 for tachyon neutrinos (B.4) 

The second version shifts the end point of the energy spectrum by m while 
the third one does not shift it. But each of the three versions has its own 
distinctive shape near the end point: The first comes in with zero slope; the 
second with infinite slope; and the third with finite slope. Earlier theoretical 
studies, as well as some experimental results, which did not survive, led to 
the hypothesis of tachyonic neutrinos. [9, 10] 

An intriguing aspect is that the Dirac equation for tachyons gives special 
prominence to the helicity of particle states. 

Also, regarding neutrinos, there are the experiments that show the effects 
of mass-mixing. The standard view is based upon the expansion of the 
wavefunction in terms of a small mass 

gi(fcx-wt) ^ e i(k(x-t)-(m 2 /2k)t) ^ (B.5) 

and then subsequent treatment of m 2 as a matrix that causes mixing with 
other neutrino states. 

For tachyons there would be only the slightest difference: 

e i(fcx-wt) ^ e i(u(x-t)+(m 2 /2w)x) ^ (B.6) 

and since these are high energy particles, u ~ k and x ~ t, it would look 
very much the same. 

One can also ask how the simplest model for cosmological expansion, with 
uniform distribution of matter, would look in the case of a tachyon-dominated 
universe. The basic equations are 

A(pa 3 ) = -3Pa 2 (B.7) 
da 



P 



const, x J E d 3 p/(e E/kT ± 1), (B.8) 



where a is the scale factor and T is the temperature, with p, P being the 
energy density and pressure as previously defined. In the case where kT » 
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mc 2 , where m is the tachyonic mass, then we have the familiar results for 
any relativistic particle: p = 3P ~ T 4 ~ a~ 4 . Now we want to see what this 
looks like in the other limit, when kT << mc 2 . Then we find that p ~ T 3 
and, if we assume Fermi statistics, then the pressure P goes as T. This leads 
to a relation between T and a as follows: 



This is a rather unusual result, predicting that the temperature goes to zero 
at some finite expansion. Of course, this very simple-minded calculation 
needs to be studied further; and it might lead to some observational com- 
parisons that help us decide whether such tachyons exist. For some related 
investigations, see [11]. 

Appendix C: Smoothed cylindrical problem 

In order to arrive at a clear solution of the field equations for the static- 
cylindrical problem studied in Section 5, we need to specify some structure for 
the source. We shall do this in a backwards manner: make up a solution for 
the metric and see what source it produces from the Einstein field equations. 
We start with a modified version of the solutions given in the earlier Section. 



as a , 



B = bs?, C = cs\ s = Vr 2 + e 2 (C.l) 
3 = a 2 /(a + 2), 7 = -2a/(a + 2). (C.2) 



We then calculate, 



_ A a{a + 4) e 2 _ C 

G °° ~ ~2B ^T2" 7' 633 " A G °° (C - 3) 



1 a 2 e 2 „ r 2 



G n G ^ ~ b Gn - (C - 4) 

Next, we integrate these densities f^r dr and find, in the limit of small a 
and small e: 



/ 



/ 



G °° " Yb a = ~ AGpl (C ' 5) 



G33 — > - Y h a = -^Gpiv 2 . (C.6) 



This verifies the results given earlier. 
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